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Abstract 

We describe an algorithm which computes components of Humbert surfaces 
in terms of Rosenhain invariants, based on Runge's method [7]. 

1 Background 

For general properties of abelian varieties we refer the reader to [2j . Denote by 
7^2 the Siegel upper half plane of degree 2, which by definition is the set 

rL2 = {T(^ Mat2x2(C) I V = T , Im (r) > 0} . 

Each r G 7^2 corresponds to a principally polarized complex abelian surface 
At with period matrix (r I2) € Mat2x4(C). Two abelian surfaces A-j- and A'^ 
are isomorphic if and only if there is a symplectic matrix M = ( ^ ^) € Sp4(Z) 
such that T = M(r) := (ar + 6)(cr + Quotienting out by this action, we 

obtain the moduli space A2 = Sp4(Z)\7-^2 of isomorphism classes of principally 
polarized abelian surfaces. It is a quasi-projective variety of dimension 3 and 
is called the Siegel modular threefold. 

The sets of abelian surfaces having the same endomorphism ring form sub- 
varieties of A2- Let j4 be a principally polarised abelian surface. Then End(j4) 
is an order in End(A) ® Q which is isomorphic to either a quartic CM field, an 
indefinite quaternion algebra, a real quadratic field or in the generic case Q. 
The irreducible components of the corresponding moduli spaces in A2 which 
have "extra endomorphisms" have dimensions 0,1,2 and are known as CM 
points, Shimura curves and Humbert surfaces respectively. 

Humbert [6] showed that for each positive discriminant A there is a unique 
irreducible Humbert surface //a in -4.2, and any matrix (t2 7-3 ) € H2 satisfying 
the equation 

kn + ^r2 - T3 = (1) 

lies on the Humbert surface of discriminant A = 4A: + £ > 0. 

The function field oi A2 is C(ji, ^'2, js) where the ji{T) are the absolute Igusa 
invariants, so there is an irreducible polynomial -f^A(iii J2) ja) whose zero set 
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is the Humbert surface of discriminant A. Unfortunately, working with Igusa 
invariants is impractical due to the enormous degrees and coefficients of the 
polynomial. One fares better by working in a finite cover of the moduli space, 
adding some level structure. Runge [7j constructed an algorithm to compute 
Humbert components in the cover r*(2,4)\7Y2 using theta functions and their 
Fourier expansions. The purpose of this paper is to apply Runge's method 
to ^2(2), the Siegel modular threefold with level-2 structure, using Rosenhain 
invariants. 

2 Level-2 structure 

Torelli's theorem implies that the map sending a curve C to its Jacobian variety 
Jac(C) defines a birational map from the moduli space of genus 2 curves denoted 
M2, to A2- Given a genus 2 curve = Y\^=iix — Ui) over the complex numbers, 
we can send three of the to 0, 1, 00 via a fractional linear transformation to 
get an isomorphic curve with a Rosenhain model: 

= x{x - l)(x - Ai)(x - A2)(a; - A3). 

The Aj are called Rosenhain invariants. 

The ordered tuple (0, 1, 00, Ai, A2, A3) determines an ordering of the Weier- 
strass points and a level 2 structure on the corresponding Jacobian, that is, 
determines a point of ^2(2). 

Let 7\/(2(2) denote the moduli space of genus 2 curves together with a full 
level 2 structure. The points of A^2(2) are given by triples (Ai,A2,A3) where 
the Aj are all distinct and different from and 1. The forgetful morphism 
A^2(2) — > A^2 is a Galois covering of degree 720 = \Sq\ where acts on the 
Weierstrass 6-tuple by permutations, followed by renormalising the first three 
coordinates to (0, 1, 00). 

As functions on A^2(2), the Rosenhain invariants generate the coordinate 
ring of 7^2(2) and hence generate the function field of ^2(2). 

3 Theta constants and Rosenhain invariants 

Let T E 7i2 and write m! = (a, b) and m" = (c, d). The classical theta constants 
(of half integral characteristic) are defined by 
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where a, b, c, d are either or 1. Write 



T?l = 


^'0000(7" 


^2 = 


^0011 (''" 


^3 = 


^'ooio(''" 


1?4 = 


^'oooi(''" 


^8 = 


^1100 (''" 


^10 = 


6'iiii(t 



These are the only theta constants we shall be using. As a function of r G A2 
there are 720 different Rosenhain invariant triples, any of which may be used. 
We use the same choice of Rosenhain triple that Gaudry uses in [3] : 

(72*^4 '^A'^IO '^2*^10 



4 Fourier series expansions 

Let us describe the Fourier expansion of even theta constants restricted to a 
Humbert surface of discriminant A = or 1 mod 4, adapted from ideas in 
Runge's paper [7j. Write A = 4k + i where i is either or 1, and k is uniquely 
determined. From equation ([T|) the Humbert surface of discriminant A can be 
defined by the set 

modulo the usual Sp4(Z) equivalence relation. Restrict Oabcd to to get 

^(r) = ^ g7ri(a;ic+a:2rf)^(2xi+a)2+fc(2x2+b)2^2(2xi+a)(2a;2+b)+^(2x2+b)2 

(xi,a;2)eZ2 

where r = e^'^*'^^/^ and q = g^'^^'^^/s^ Unfortunately, q has negative exponents 
which computationally makes it difficult to work with this expansion. To over- 
come this difficulty, make the invertible substitution r = pq to produce the 
expansion 

^_-^^-^xic+X2dp(2x-i+a)'^+k{2x2+b)^^{2xi+a+2x2+b)'^+{k+e-l)(2x2+bf 



(xi,x2)e: 



7.2 



which is more computationally friendly, being a power series with integer coef- 
ficients. Call the above expansion the Fourier expansion of Oabcd restricted to 
Ha. 

Addition and muliptication of restricted Fourier expansions are just the 
usual addition and muliptication operations in Z[[p, g]]. To compute the ex- 
pansions of Rosenhain invariants we need to know how to invert elements of 
Z[[p, g]] where possible. It is well known fact about power series rings that if 
f{p,q) is in Z[[p, g]] with /(0,0) 7^ , then f{p,q) is a unit with inverse given 
by the geometric series 



3 



An implementation on a computer uses truncated Fourier expansions, where 
arithmetic is done in g^) for some positive N. It is easy to see 

that the geometric ratio has zero constant term, in particular (1 — ///(0,0))'^ € 
(p^, q^) for k > N so the above formula converges to the truncated expansion 
q£ j-i £qj, g^j^y g]2osen precision. 

From the expansions we observe that 'i^i, ??2i ''^a, ''?4 have constant term 1, 
hence are invertible, but i^s = 2p^^^q^^^~^ + . . . and t?io = —2p^^^q^^^~^ + . . . 
have zero constant term. Fortunately one can show that i?8,i?io are in the 
ideal {p^^^q^^^~^)'L[[p,q\] hence by cancelling out the p^+^^q^+^-'^ factors, the 
quotient t^s/^io makes sense in g]]. Thus we are able to compute the 
Rosenhain invariants Ai,A2,A3 as Fourier expansions restricted to a Humbert 
surface. 



5 The algorithm 



We describe an algorithm to find the equation of all irreducible components of 
T^A representable in terms of Rosenhain invariants. 

Much arithmetic-geometric information is known about Humbert surfaces, 
and more generally Hilbert modular surfaces (see [5], [8]). All Humbert compo- 
nents are hypersurfaces in ^2(2) of the same degree. Of use for us is the degree 
m(A) of the defining (irreducible) polynomial -Fa,? of a Humbert component as 
well as the number of components of H^. The number of Humbert components 
in the Satake compactification A^i^) of ^2(2) is (see [Ij) 



m(A) 



10 if A = 1 mod 8 
15 if A = mod 4 . 
6 if A = 5 mod 8 



The degree of an irreducible polynomial ■ defining a Humbert component 
in ^2(2) is given by a recursive formula 

aA = ^m(A/x2)deg(F*^/^2) J 

where a a is the coefficient of a certain modular form of weight 5/2 for the group 
ro(4), which fortunately has a more elementary description due to a formula 
of Siegel, 



OA - 24 ^ cTi 




x'^\ J 12A — 2 if A is a square 
otherwise 



This provides an upper bound on the degree of the polynomials i^A,i(ei, 62, 63). 
From computational evidence it appears degFA = degF^ for nonsquare dis- 
criminants A and that degF„2 = (1 — ^) degF*2 for all n. 

The algorithm is very simple. We have 61,62,63 represented as truncated 
power series. We know the degree of the relation we are searching for. To 
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A 


1 


4 


5 


8 


9 


12 


13 


16 


17 


20 


21 


24 




1 


4 


8 


8 


24 


16 


40 


32 


48 


32 


80 


48 



Table 1: Table of degrees 

find an algebraic relation of degree d, compute all monomials in ei , 62 , 63 of 
degree atmost d and use linear algebra to find linear dependencies between the 
monomials. Once one component has been determined, the others can easily 
be found by looking at the Rosenhain S'e-orbit of a component. 

The fixed groups of the Humbert components in this model can be com- 
puted. As we know, Sq acts on the Rosenhain invariants via the natural action 
on (0, 1, 00, ei, 62, 63). Let /ia be the Humbert component computed using the 
above algorithm. The fixed group of /ia for even discriminant splits into two 
cases, 



Fixsg(/i4fe) 



G if A: is odd 

g^^Gg if k is even 

where G C is a group of order 48 generated by three elements 

(0,61,63,00,62,1), (61,62) and (1,61,63,62); 

the conjugating element \s g = (1, oo)(6i, 62, 63). Ignoring discriminant 1 which 
is a special case, the fixed group of A = 1 (mod 8) is a group of order 72 
generated by 

(0, 6i)(l,62)(oo,63), (1,00), (61,62) and (62,63). 

For A = 5 (mod 8) the fixed group is a group of order 120 generated by 

(0,6i)(l,62)(oo,63), (1,63,62,61,00) and (00,61,63,62). 

By making use of some of the simpler fixed group symmetries, we can reduce 
the size of the linear algebra computation. For example, the discriminant 12 
component /112 satisfies /ii2(e2, ei, 63) = /ii2(ei, 62, 63) which means we only 
need roughly half the number of evaluated power series since efe'^e'^ and eXe'^e'^ 
have the same coefficient. 

The runtime of the algorithm is greatly affected by the 0(('^3'^)) = 0{d?) 
monomials that need to be evaluated. The linear algebra solution requires 
finding the kernel of a matrix with 0{{^~^^^) rows and in the order of (A^/4)^ 
columns where N is the precision of the power series, which gives a runtime 
cost is 0{d^N'^). To have any chance of finding a unique relation, the number 
of monomials must be less than the precision used, so that the runtime is at 
least of order 0{d^). 

From the table it is evident that the degree increases with the discriminant, 
so as it stands this algorithm can only find equations with small degrees. Be- 
sides discriminant 21, we managed to produce Humbert components for all the 
discriminants listed in the above table. This extends the equations found in 
the literature ([6],[1]) which go up to discriminant 8. See the appendix for the 
equation of a discriminant 12 Humbert component we found. 
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Appendix: Equation for discriminant 12 

= e^el - 4e|e| + 6e|e^ - 4e|e| + efel - ieie^el - 16eiefe| + 
40eieie|-16eieie|-4eieie| + 160eie|e|-160eie|e|-160eie|e| + 
leOeiefe^ - 132eie|e| - 272eie|e| + 808eie|e| - 272eie|e| - 
132eieie^ + 384eie|e| - 384eie^e| - 384eie^e| + 384eie|e^ - 
256eie^e|+512eie^e|-256eie^e|+6efe|e|+40efe^e| + 164efe|e^ + 
40efeie^+6efe^e|-160efeie|-352efeie|-352efeie^-160efeie^- 
272efete| + 1344efe|e| - 608efe^e| + 1344e?e^e^ - 272efe^e^ + 
384efe|e| - 416efe|e| - 480efe|e| - 480e?e|e| - 416efe|e^ + 
384efe|e| - 762e2e|e| + 1064e2e|e| - 348efe|ei + 1064efe|e| - 
762efe^e^ + 384efe^e| - 384efe^e| - 384efe5e| + 384efe^ei - 
4ef 626^ - 16ef e2e| + 40ef 626^ - 16ef 626^ - 4e?e2e| - 160ef e|e| - 
352efeie| - 352efe|e| - 160efe^e^ + 808efe|e| - 608efe|e| + 
3696efeie| - 608efefe^ + 808efe|e_^ - 384efe|e| - 480efe|e| - 
2208efe|e| - 2208efe^e| - 480efe|e^ - 384efe|e5 - 256efe|e3 + 
1064efe|e| - 608efe|e| + 3696e?e|e;^ - 608e?e|e| + 1064efe|e^ - 
256efe^e^ + 384efe|e3 - 416efe|e| - 480efe|e| - 480efe^e| - 
416efe^e| + 384efe|e| - 132efe^e3 - 272efe|e| + 808efe^e| - 
272efe^e| - 132efe^e| + efej - 4efe| + 6efe^ - Aefel + efe| + 
160efe2e| - 160efe2e| - leOefeae^ + 160efe2e^ - 272efeiei + 
13Uefelej - 608efe2e| + 1344efe|e| - 272efe2e| - 384efefe| - 
mefelel'- 2208efe|e| - 2208efe|e| - 480efeie^ - 384efe|e5 + 
512e^e|e3 - 348efe|e| + 3696ef e|e| + U96efejej + 3696ef e|e| - 
SASefejel, + 512efe|e^ - 384ete^e3 - 480efe|ei - 2208efe|e| - 
2208efe|e| - 480e|e|e| - 384efe|e| - 272efe^e3 + 1344efe|e| - 
608efe^ei + 1344efe|e| - 272efe^e| + WOefeles - 160efe^e| - 
160efe^e| + 160efe^e| + e|ef - 4efefe3 + 6efe|e| - 4efe|e| + 
efe|e|-132efe2ei-272efe2e|+808efe2e|-272efe2e^-132efe2e|+ 
384efe^e| - 416efeie| - 480efeie| - 480efe^e| - 416efe|e^ + 
384efe^e^ - 256efeie3 + 1064efe|e| - 608efe|e| + 3696efefe| - 
608efe^e| + 1064efeie^ - 256efe^e| - 384efe|e3 - 480efe|e| - 
2208efe|e3 - 2208efe^e| - 480efe|e| - 384efe^e^ + 808efe|e3 - 
608efe|e|"+ 3696efe|ei - 608efe|e| + 808efe|e| - 160efe^e3 - 
352ef e^el - 352ef efef - 160ef e|e| - 4efe^ - 16ef 6^63 +40ef e^e^ - 
16efe|ei-4efe^e|+384efe2ei-384e5e2e|-384e'i*e2e|+384e5e2e^- 
762efeie| + 1064efeie| - 348efe|e| + 1064efe|e| - 762e^e^e| + 
384e?eie3 - M^e\4,e\ - ^me\e\el - ^me\e\e\ - 416efe|e| + 
384efe|e^ - 272efe|e3 + 1344efe|e| - 608efe^e| + 1344efe|e| - 
272e?e|e|-160e^e|e3-352e^e|e|-352e?e|ei-160efe|e|+6e'i*e| + 
40e?e|e3 + WAefe^ej+AOefe^el, + 6e'i'e|e| - 256e[e2e| + 512e[ e2e| - 
256e[e2e| + 3846^6161 - 384e[eie| - 384ejeie| + 384ejeie| - 
132eje3e3 - 272eje|e| + 808eje3e| - 272e3'e3e| - 132e^efe| + 
160eje^e3 - 160eje|e| - 160eje|e| + 160eje|e| -4ej'e| - 16eje|e3 + 
40eje|e| - 16eJ e|e| - 4e[e|e| + e^e^ - 4ef e|e3 + 6ef e|e| - 4e?e|e| + 
ef 6263 
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